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I consider the qnantity 5{mim2m^) = ~ (-^9192 + -^9293 + ^qiqa)/^^ where 

the M’s represent the gronnd state spin-averaged hadron masses with the indicated qnark 
content and the m’s the corresponding constitnent quark masses. I assume a logarith¬ 
mic interquark potential, the validity of a nonrelativistic approach, and various stan¬ 
dard potential model inputs. Simple scaling arguments then imply that the quantity 
R{x) = d{mmms)/d{mQmorno) depends only on the ratio x — m/ms, and is independent of 
mo as well as any parameters appearing in the potential. A simple and accurate analytic de¬ 
termination of 5{mmm^), and hence R{x), is given using the 1/D expansion where D is the 
number of spatial dimensions. When applicable, this estimate of R{x) compares very well 
to experiment — even for hadrons containing light quarks. A prediction of the above result 
which is likely to be tested in the near future is M-£*/2 + {Ma^+Mj:,^)/4 = 5774±4 MeV/c^. 
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In 1983, Nussinov suggested that 

^ i^qiq2 "I" -^5293 ^qiqa)/^ (1) 

is a rigorous inequality in QCD. Here Mq^q^q^ represents the mass of any spin state of the 
ground state baryon with the indicated quark content, and Mq^q. the mass of the ground 
state meson in the spin state imposed by the chosen baryon.0 He was able to demonstrate 
this assuming a separation of the three-body Hamiltonian for the quarks in the baryon 
into a sum of two-body pieces — a result which holds at both weak and strong coupling 
in QCD. Although this does not constitute a complete proof, the inequality (|^) is indeed 
satished by the observed ground state baryons and mesons. The arguments in |]^ apply 
not only to specihc spin states as noted above, but also to spin-averaged states . That is, 
we can take the M’s in (|l|) to represent the spin-averaged ground state baryon and meson 
masses. This view of the Nussinov inequality, which we adopt below, is also supported by 
experiment. 

Is it possible to further understanding the relationship between ground state baryon 
and meson masses? More precisely, if we write 

92 93 ~ (Mgiqrj -|- Mq^q^ Mq^q^')/‘I -|- 6 (17117712^11^) , (2) 

can we say anything about the positive quantity 5(mim2ms)l How does 5 depend on the 
quark masses mi, m2 and m3? Does anything special happen when one or more of these 
masses get very large, or when one or more are equal? To begin to answer these questions 
clearly requires more specihc dynamical input than that used in arguing for the inequal¬ 
ity (|l|). A hrst principles analytic computation from QCD is clearly beyond our current 
theoretical resources. Numerical results from lattice QCD have not yet reached the level 
of accuracy required to address these questions seriously. So, in order to proceed at all, a 
more phenomenological approach is called for. For instance, we can attempt to study 5 in 
the context of potential models. That is, we can postulate a specihc form for the interac¬ 
tion between the quarks in mesons and baryons, and then solve the appropriate two and 

^ This inequality was also discussed independently in the context of nonrelativistic potential 
models by Richard [^. The special case of equal quark masses was hrst treated by Ader, Richard 
and Taxil |Q. Similar inequalities were also derived at this time by Weingarten |Q and Witten 
||. Various “improved” versions of (|) have also been discussed over the years. See, for example, 
iP and references therein. 
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three body problems to obtain the desired spectrum. Of course, there are necessarily free 
parameters in this approach which must be ht to a subset of the experimental data. This 
is because we do not know the exact form of the interquark potential, nor do we know the 
exact constituent quark masses that should appear in the above Hamiltonians. Moreover, 
three body equations are notoriously difficult to solve or approximate analytically, so that 
one must in general resort to sophisticated numerical methods. So even this very natural 
idea of potential models, although having been applied very successfully for over 30 years, 
is found to lack a certain simplicity and elegance. 

We will hnd, however, that a simple, accurate, and completely analytic approximation 
to 6 can be obtained if one assumes a logarithmic interquark potential and nonrelativis- 
tic dynamics governed by the Schrodinger equation — at least when two of the quark 
masses are equal. The log potential approach to the hadron spectrum has been around 
for two decades [0. It was originally motivated by the observation that the '0(2S')-'0(1 S') 
and T(2S')-T(1S') mass splittings are approximately equal, as this splitting is quark mass 
independent in the log potential. After many years, and the discovery of many new states, 
hts to the mass spectrum of mesons containing b, c and s quarks still support the quasi- 
logarithmic nature of the interquark potential ^j. There is also evidence supporting this 
picture in the low-lying baryon spectrum, although computations are more difficult here 
and experimental results more sparse [Q0. Of course, it is clear from QCD that this 
description must eventually break down at very short and very large distance scales. How¬ 
ever, more realistic QCD-inspired potentials — such as the famous “Coulomb-plus-linear” 
potential [ITO — are all quasi-logarithmic over the distance scales relevant for the low-lying 


spectrum of the observed heavy hadrons [^.H In short, the log potential and its relatives 
have been very successful. 

The approach used here is to approximate 5 to leading order in the 1/D expansion, 
where D is the number of spatial dimensions. That is, the relevant two and three body 
Hamiltonians are hrst generalized from 3 to T> spatial dimensions, and then we assume that 
D is large. It is well known that a systematic 1/D expansion can be developed [1^] |T4 


There is also strong evidence that the leading order term is already very accurate for quasi- 
logarithmic potentials, in both the two and three body cases, even for D = 3. The two 


^ Other physical quantities such as production and decay rates, which are sensitive to the wave 
function at the origin, depend on physics at shorter distance scales. Here there are substantial 
differences between the predictions of quasi-logarithmic potentials and more realistic ones, even 
though they have the same “shape” near the RMS radii of the states. See, for example, [ 1 ^|. 
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body results in the l/il> expansion also have the nice property that they are completely 
analytic for all power-law potentials V{r) = Ar^. Since £n{r) = hmb^o(?'^ ~ 1)/^: the 
results are also analytic for the log potential. In the three body case this property no 
longer holds for general power-law potentials except in the case of three equal masses. For 
instance, when only two masses are equal, the leading order large-result requires the 
numerical solution of a transcendental algebraic equation [^. However, for the log limit I 
will show that a completely analytic solution can be found. Moreover, when put together 
with the two body result, this leads to an extremely simple and accurate form for S. (For 
three unequal masses the situation is more complicated |]l^ , and will not be treated here. 
Related applications of large-results to the ground state meson and baryon spectra can 
be found in [P^[|TB| .) 

Lets begin with the specihc form of the interquark forces. I will assume that the po¬ 
tential between a quark and an antiquark in a color singlet meson is Vg^g^{r) = Ainirjr^, 
where r is the interquark separation and H, ro > 0 are parameters. I make the further stan¬ 
dard assumption that the quark-quark potential in a baryon (that is, two quarks in an 
overall color antitriplet) is 1/2 of the color singlet quark-antiquark potential at equal sepa¬ 
ration: Vg^g^ (r) = Vg^g^ if) j ' (l iguoro the effects of the direct three body interactions that 
arise in QCD, as these have been shown on many occasions to be negligible for the ques¬ 
tions that we are considering. See, for example, iOO I also assume that the dynamics 
are governed by the appropriate nonrelativistic Schrodinger equation. If 1^2'?/’ = Errumj'^ 
and Hs(f) = cire the two and three body Schrodinger equations (for the ground 

state) respectively, then we have 5{mim2m’i) = Em^rn 2 m 3 - {Em-,m 2 +-£^mim 3 )/ 2 . 

Throughout, the effects of spin are neglected and the results are interpreted as representing 
the masses and binding energies of spin-averaged meson and baryon states. The simple 
rescaling TiVi/^/rriiA of the quark position vectors in H 2 and shows that 


E. 


miTTij 


E, 


mi m2 777,3 


{A/2)[in{h'^/nArl) -F C] 

= {3A/4)[in{h‘^/miArl) + f{w,x)], 


(3) 


where w = mi/m 2 , x = mi/m^, and fx = mimj /(m^ -1- mj) is the two body reduced mass. 
Both the constant C and the function / are independent of h, A and tq. These scaling 
relations further imply that 

5{mim2m^) = Ag{w,x), (4) 
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where again g is independent of any dimensionful parameter in the problem other than the 
constitnent qnark masses, and even this dependence is only throngh the two dimensionless 
ratios w and x. 

Note that an immediate conseqnence of (^) is that in the case of eqnal qnark masses 
mi = m 2 = m 3 = m, 5 is independent of m. Assnming that the constitnent masses of the u 
and d qnarks are identical and eqnal to some valne m^, where here and below the symbol 
n denotes either u or d, we may then determine 5{mnmnTnn) from experiment. More 
precisely, 5{mnmnmn) — Mnnn — ^{Mnn)/2. Assnming fnrther the standard techniqne 
for constrncting the spin-averaged masses, we have Mnnn = [Mn + M/^)l2 and Mnn = 
(Myr + 3Mp)/4. Pntting this all together with the measnred masses! of the A, tt and 
p states gives dimn'rrin'n^n) = 171 ± 2 MeV/c^. Thus, we should expect 5{msmsms) to 
also be in this range, where mg is the mass of the strange quark. Unfortunately, this 
cannot be tested directly since the spin 1/2 ground state baryon with three s quarks does 
not exist. There is only the spin 3/2 0“. (This is because the spin wave function of 
the identical s quarks must be totally symmetric in the ground state.) There is a similar 
problem in the meson sector where only the spin 1 (j) can be interpreted as a nearly pure 
ss ground state. (The spin 0 combination mixes strongly with uu and dd states.) Thus, 
some additional theoretical input is needed to obtain the spin-averaged baryon and meson 
masses. However, any reasonable model of the spin-spin interactions between quarks, when 
applied to the 0 “ and (j) in order to reconstruct the spin-averaged states, will yield a value 
of 5{msmsms) around the desired range. A similar test of d{mmm) for charm and bottom 
quarks cannot be made since the heavy ccc and bbb baryons have not yet been observed. 

It is interesting to note that a numerical evaluation of 6 in the case of equal quark 
masses yields 5{mmm) ~ (0.225) A |p. (That is, g{l, 1) ~ 0.225.) If we then use this result 
to “ht” A to the observed 5 for u and d quarks, we obtain A = 760 ± 10 MeV. This is not 
too far from the value A = 733 MeV 0 obtained from a ht of the log potential spectrum 
to meson states containing only heavy quarks! Not only do we see a striking consistency 
between the meson and baryon sectors, but we also gain some faith in the application of 
our nonrelativistic model to hadrons containing “light” quarks. This should not be too 
surprising since nonrelativistic models have often been successful on such semirelativistic 
quark systems in the past. The standard philosophy regarding this is that many relativistic 


^ Unless otherwise stated, all experimental results are taken from [1^|. I also average over the 
electromagnetic splittings within a given isospin multiplet. 
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corrections can be incorporated into a nonrelativistic model by a redefinition of the con- 
stitnent quark masses and simple modifications of the potential such as the addition of an 
overall constant, at least for low-lying states H- (Note that both of these modifications 
do not ahect the result for 5{rnmm) in the log potential.) Moreover, in our case there 
may be some cancellations of additional relativistic corrections occurring between the two 
body and three body contributions to 5. I will maintain this viewpoint in what follows 
and continue to apply the results of the above model to systems containing light and/or 
heavy quarks. 

What about the case when only two of the quark masses are equal? Here, 5{mmm^) = 
Ag{l,x), where x = mlm^. Instead of resorting to a time-consuming numerical determi¬ 
nation of as a function of x, lets consider 5 to leading order in the 1/D expansion. To 
this end, let and (j) — be the generalizations of the two 

and three body Schrodinger equations to an arbitrary number of spatial dimensions. That 
is, we keep the interquark potentials the same and simply replace each three-dimensional 
Laplacian appearing in the kinetic energy terms by its T>-dimensional counterpart. We 
can still write EmXij Em^rn 2 ra 2 . as in (^, only now the constant C and the function 
f{w, x) depend on D. These binding energies can be obtained in the leading order of the 
1/D expansion by taking the appropriate limit of the corresponding results for power-law 
potentials. For the two body case, this has already been done |2^. The result is 


= (.4/2)|fa(efi"l)V4MrS) +0(1/1))]. 


(5) 


The T> —> cx) limit is a classical limit (distinct from the —> 0 limit), and the above result 
can be interpreted as I4ff(-Ro) where V/ff is the large-T> ehective potential. Here Rq is the 
interquark distance at which this ehective potential is minimized. In other words, Rq is 
the “size” of the meson. For the above log potential we have Rq = \J?i^/AgA. 

In the three body case, the leading order large-T> result for power-law potentials is 
not completely analytic, as noted previously. However, the log limit can still be recovered 
in closed form when two of the masses are equal. A careful analysis yields 


dW, = (3>l/4){fa|eR"D"(l + 2x){2a)^lg2raArl(2 - o)"] + 0(1/D)}, (6) 

where x — m/vri'i and a = 8/[^(2x -|- l)(2x -|- 25) -t- 2x -t- 5]. This corresponds to a 
minimum energy configuration of the large-T> ehective potential in which the distance 
between qj, and either of the other two quarks is W h^(l-|-2x)/2mA(2 — a)^, and the 
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angle Q between the two associated displacement vectors satisfies cos(6*) = 1 — a. Af¬ 
ter some messy bnt straightforward algebra, one obtains the following expression for 
5{mmms) = E^Xm^ - EXXI 3 “ eXXII‘2: 

5 {mmm 2 ,) = (A/4) [Mn{\/{2x + l)(2x -|- 25) -|- 2x -|- 5) -|- 3£n{2x -t- 1) 

/- 

{2x + l){2x + 2h) + 2a: + 1) - 2£n{x + 1) + 0{l/D)]. 


Note that 5 {mmm^) is independent oi D D ^ 00 . When a: = 1, that is m3 = m, 
we have S{mmm) = (3A/4)£n(4/3) ~ (0.216)A, which compares well with the n um erical 
resnlt of (0.225)A. (It is interesting to note that the large-result for 5 {mmm) is exactly 
the same as that obtained in H by seemingly independent means.) It seems reasonable 
to believe that at least some of this accuracy is due to a partial cancellation of the 1/D 
corrections between the two and three body cases. We can also easily extract the results 
for X = 0 and x ^ 00 from (^. We have 5(x = 0) = (A/4)£n(3125/1458) ~ (0.191)A and 
d{x 00) = (A/4)£n(2) ~ (0.173)A. Scaling arguments imply that the latter result is 
actually exact! Unfortunately, in this limit we do not expect the log potential to provide 
an accurate description of quark dynamics. This is because we can think of the x ^ 00 
limit as either m getting very large or m3 getting very small. In the former case, the two 
very heavy quarks of mass m in the baryon are clearly in the perturbative “Coulombic” and 
not the logarithmic regime, while in the latter case a nonrelativistic treatment of the very 
light quark is not justified. (Similar problems do not occur in the x —0 limit viewed as m 3 
getting very large.) As mentioned earlier, however, we do believe that the nonrelativistic 
log potential is applicable to systems containing u,d,s,c and b quarks. Here, x can get as 
large as mb/mn — 15. 

It is often a bit more useful to consider the ratio R{x) = (j(mmm 3 )/fi(momomo). This 
quantity has the advantage of being independent of the parameters appearing in the log 
potential. A plot of R{x) in the large-17 limit is shown in Fig. 1. It possesses a single 
critical point (a global maximum) in the range 0 < x < cx). This occurs at x = 1 where 
i?(l) = 1. (Note that the large-17 result for i?(x) is forced to be exact at x = 1, so 
that one might expect it to be somewhat more accurate in the region of physical interest 
than the corresponding result for fi, which already seems to be within a few percent of 
the exact number for all x. More on this below. Of course, R{x) is no longer exact as 
X —> cxo. But this is a trade-off that we will gladly accept since, as noted above, our model 
is not applicable at very large x.) The asymptotic values 72(0) ~ 0.883 and R{oo) ~ 0.803 
represent completely parameter independent relations for the log potential, at least in the 
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large-D limit. Note also that R{x) only changes by about 20% over an inhnite range of 
X values, and only by about 10% in the physical region mn/mb < x < mb/rrin- Higher 
order 1/D corrections to the two and three body energies, and hence to 5 and R, can in 
principle be computed. Although the resulting formulas will be more complicated, it is 
nice to know that there is a systematic analytic procedure for improving our results. For 
instance, using the results of one obtains 

5{mmm) = (3A/4)[£n(4/3) + (2^3 - ^2 - 2)/H + 0{l/D‘^)] 

( 8 ) 

= [3£n(4/3)A/4]{l + a/D + 0{1/D‘^)}. 


For D = 3 this yields d{mmm) ~ (0.228) A which is 99% accurate, to be compared to the 
96% accuracy of the leading term alone. Of course, the above 1/D correction has no effect 
on i?(l) which is already exact at leading order. However, for a; % 1 our result for R{x) 
will have 1/D corrections. As an example, consider i2(0). A lengthy computation using 
the results of yields 


S{x = 0) = (A/4) [£n(3125/1458) + 6(v^ + ^17 + 2^6 - 5^2 - 5)/5D + 0(1/D^)] 

= [£n(3125/1458)A/4]{l + P/D + 0{1/D^)}. 

(9) 

For D = 3 this gives d{x = 0) ~ (0.202)A, compared to the leading order result of 
approximately (0.191)A. Putting (H) together with (H) we obtain 


R{0) = [£n(3125/1458)/3£n(4/3)]{l + (/? - a)/D + 0{1/D^)}. (10) 


For D = 3 we have {P — a)/D ~ 0.0016, so that the change in i?(0) due to 1/D corrections 
is about 1/6 of 1%. (i?(0) changes from about 0.8833 to approximately 0.8848.) This is 
quite remarkable convergence. So although not exact as for i?(l), the leading order result 
for i2(0) should be extremely accurate. I believe the above examples to be illustrations of a 
quite general pattern: the leading order large-H result for 5{mmm^) is accurate to within 
a few percent, while the corresponding result for R{x) is accurate to within a fraction 
of a percent — at least over a large range of x. To be safe, one can compute the 1/D 
corrections to R{x) for a general x using the results of [^. However, it seems clear that the 
leading order results are accurate enough in general that the added complexity introduced 
by adding higher order corrections is not justihed. 
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Figure 1: R{x) in the large-limit. 

How do the above results compare with observation? Lets hrst consider x = mnjms 
and X = msjrnn (we have already discussed the situation for x = 1). In the former case 
we have 

5{mnmnms) = (Ma My)/A + My*/2 

- (M^ + 3M,)/8 - {Mk + 3Mk*)/A (11) 

= 171±2 MeV/c4 

Recalling that 5(mnm„mn) = 171 ± 2 MeV/c^ as well, this leads to the experimental 
value RexpiTTT'n/TTT's) = 1.00±0.02. However, we cannot obtain an unambiguous theoretical 
prediction for R{mn/ms) from (|^) since and ms have not been determined in our model. 
One can use, among other things, the ratio of the K-K* and n-p mass splittings as a naive 
estimate of This gives m^/mg ~ 0.63. (For comparison, the ratio of the E-E* to 

the N-A splitting yields mnjms — 0.65.) But luckily, R does not change too much within 
the entire range of “reasonable” values of mnjms- For 0.5 < mnjms < 0.8, R{mnjms) 
stays between 0.985 and 1.00, in good agreement with the experimental number. Using the 
corresponding range 1.25 < mg/m^ < 2.0, the log model prediction for i?(mg/m„) also 
varies between 0.985 and 1.00. Unfortunately, we cannot obtain a purely experimental 
number to compare this with since certain states needed to hnd the required spin-averaged 



masses Mgs and Mgsn do not exist. However, estimating the masses of these missing states 
from other spin-splittings of hadrons containing u, d and s qnarks yields an “experimental” 
valne of R{ms/mn) in the desired range. 

We now tnrn to hadrons containing one heavy qnark — that is, one c or b qnark. 
Taking the ratio of the D-D* and n-p mass splittings yields the estimate mn/rric — 0.22. 
(Using the recent measurement M-£* = 2519 ±2 MeV/c^ made by the CLEO collaboration 
||24|| , we can also estimate rrin/mc from the ratio of the — E* and N — A splittings. 
This again gives mn/rric — 0.22.) From we then obtain the prediction R ~ 0.95. As 
before, the deviation from this value is not too large over the entire range of reasonable 
values of mn/rric- More precisely, for 0.15 < rrin/mc < 0.3, the value of R varies from 
approximately 0.935 to 0.965. The experimental formula for 5{mn'mnmc) can be obtained 
from (O) by replacing K and K* by D and D*, as well as substituting Ac, Ec and 
E* for A, E and E*. This yields 5{mnmnmc) = 165 ± 2 MeV/c^, which corresponds to 
i?exp(™n/wc) = 0.965 ±0.02, once more in good agreement with the theoretical prediction. 

Similarly, the ratio of the B-B* and 7r-p splittings yields mn/vrih — 0.07, which when 
plugged into gives R ~ 0.91. Moreover, for the range 0.05 < mn/mi, < 0.08, the value 
of R stays between 0.905 and 0.915. The experimental value can again be obtained from 
(pTl), this time by substituting B and B* for K and AT*, and replacing A, E and E* by Af,, 
Ef, and E^. But since the experimental status of some of the bottom baryons is uncertain 
(see below), I instead use this formula in conjunction with the above result for R to obtain a 
rather precise prediction for the bnn center of gravity — namely, Mi,nn = 5774±4 MeV/c^. 
To what extent can this result be checked? In other words, what is our current state 
of knowledge about the bnn baryons? As of a year ago, the A}, had been seen but its 
mass was not known very accurately, the world average being 5641 ± 50 MeV/c^. Since 
then, there have been three additional measurements. The most accurate is that of the 
CDF collaboration |^. They found a Af, mass of 5621 ± 5 MeV/c^. The other two 
results are from the ALEPH and DELPHI collaborations. The ALEPH number is 
5614 ± 21 MeV/c^, while DELPHI found 5668 ± 18 MeV/c^. Taking the weighted average 
of these new results and the previous world average gives a new world average for the Af, 
mass of 5624 ± 5 MeV/c^, which we will use below. Recently, the DELPHI collaboration 
has also reported the hrst evidence for the Ef, and E^ baryons [^. Their preliminary mass 
values^ are — Ma,, = 173±14MeV/c^ andM^*—M a,, = 229±9MeV/c^. Noting that 


I thank C. Kreuter for informing me of the new (larger) error for — Ma,,. 
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Mbnn = + {My,* — )/2 + {My,, — M ^^) /4, we can plug in the above measurements 

and obtain the experimental value M^^ = 5782 ± 8 Mev/c^. Though consistent with the 
above prediction, it would clearly be desirable to have a more precise measurement of all 
the bnn baryon masses in order to make a better comparison. At this point it is also worth 
noting one strange feature of the measured values of the "Ef, and masses — namely, their 
difference.! We obtain from the DELPHI numbers that Me* — Mej, = 56 ± 17 MeV/c^. 
An estimate of this difference from the spin-splitting pattern in the meson sector would 
be around 20 MeV/c^. (By contrast, recall the amazing consistency of the meson and 
baryon hyperhne splittings in the snn and cnn cases.) Various other simple estimates 
yield a similar number Given this discrepancy, I would not be surprised to see 

these preliminary mass values change somewhat in the near future and come into better 
agreement with both the mesonic spin-splittings and the predictions of the log model. 

It is worth stressing again that there is every reason to believe the above predictions to 
be extremely trustworthy. As noted earlier, the nonrelativistic log model should provide an 
accurate description of the ground state hadrons for some values of the parameters, even for 
light quarks. However, all of the parameter dependence has cancelled out in R{x) except 
for a relatively weak dependence on the quark mass ratio x, and we have incorporated 
our ignorance about these ratios into the overall uncertainty. Finally, 1/D corrections also 
seem to be completely under control — the additional uncertainty in our results due to the 
neglect of such corrections is swamped by the quark mass uncertainties just mentioned. The 
appeal of the above approach is in its simplicity, accuracy and parameter-independence, 
and not in any rigorous connection to more fundamental ideas. 
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